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I. INTRODUCTION 

Recent several experiments require the serious modifications of the standard model (SM) 
in spite of the great success. For an example, the SM Higgs boson search in the diphoton 
mode h — > 77 at Large Hadron Collider (LHC) is shown that its signal strength is 1.65 ±0.24 
at ATLAS jl, 2] and 1.6±0.4 at CMS For another example, the existence of non-baryonic 
dark matters (DMs), which cannot be included in SM, dominates about 23% from the CMB 
observation by WMAP Q. The fact is strondv supported by the cosmological observations 
such as the rotation curves of the galaxy [5[ and the gravitational lensing [6[ in our uni- 
verse. In recent years, direct detection experiments of DM; XENON100 CRESSTII js], 
CoGeNT |9] and DAMA 10 1 , show the scattering events with nuclei. XENON100 has not 
shown a result of DM signal but shown an upper bound with the minimal bound around 100 
GeV. On the other hand, CoGeNT, DAMA and CRESSTII have reported the observations 
which can be interpreted as DM signals that favor a light DM with several GeV mass and 
rather large cross section. As far as we consider these experiments, the mass scale of DM 
should be 0(1-100) GeV. 

In order to explain the excess in the h — > 77 channel, a modified diphoton event has 
been discussed in a lot of paper. If a model contains charged new particles which couple 
to the SM-like Higgs boson such as charged Higgs bosons, the decay rate of h — > 77 can 
be enhanced due to loop effects of these charged particles. However, in a model with only 
one pair of singly-charged scalar bosons such as two Higgs doublet models, it is difficult to 
predict around 60% enhancement of the decay rate unless the mass of the charged scalar 
bosons is taken to be smaller than about 100 GeV 1 . The minimal Higgs triplet model (HTM) 



12j, introducing the isospin triplet scalar 



motivated from the type II seesaw mechanism 
field A, can easily explain the diphoton anomaly. Moreover, if A can be inert scalar, then 
its neutral component can be a promising DM candidate that can strongly correlate with 
neutrinos. This is like a radiative seesaw models 13l442l|. Hence the theory can be realized 



at TeV scale, so be well-tested at current experiments like LHC. 

In this paper, we propose a two-loop induced neutrino model with gauged B — L symme- 



1 If the decay rate of the Higgs to bb mode is sufficiently suppressed compared to the SM value, the branching 
ratio of the Higgs to diphoton mode can be enhanced without changing its decay rate. In that case, the 
branching fraction of the Higgs to WW* and ZZ* modes are also enhanced . 
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try that is an extension of the HTM. In the bosonic sector, three scalar fields are introduced 
in addition to the SM particles; A and t] are an SU(2)i triplet and doublet fields, respec- 
tively which do not have vacuum expectation values (VEVs), and x is SU(2)l singlet which 
acquires a VEV after the spontaneous B — L symmetry breaking 2 . In the fermionic sector, 
we introduce an exotic vector-like lepton with SU(2)l doublet 44J, and three right-handed 
neutrino with SU(2)l singlet, both of which can contribute to the radiative neutrino mass. 
Due to the abundant extra fields, we have several DM candidates. Also we can discuss the 
testability of the Higgs sector, especially, the doubly-charged scalar boson, in which we could 
have a discrimination to the HTM, since it couples to the exotic lepton. As a notice, such a 
complicated model can be realized within non-Abelian symmetries. So we briefly show how 
to realize our model in the appendix. 

This paper is organized as follows. In Sec. II, we show our model building including 
the Higgs potential, stationary condition, neutrino mass, and lepton flavor violation (LFV). 
In Sec. Ill, we analyze DM phenomenologies. In Sec. IV, we analyze Higgs phenomenology 
including electroweak precision observables, signatures of the doubly-charged scalar boson at 
LHC We summarize and conclude in Sec. V. In appendices, some result of detail calculations 
are write down; gauge boson two point functions, decay rate of the doubly-charged scalar 
boson, and the assignments for each particles in non-Abelian symmetries. 



II. THE RADIATIVE SEESAW MODEL 



A. Model setup 
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TABLE I: The particle contents and the charges for fermions. V and L' c are exotic leptons. 



2 There exists a B — L gauge boson, however we neglect through our analysis because it can decouple to 
the other particles from the LEP II 43 1. 
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TABLE II: The particle contents and the charges for bosons. 



We propose a two-loop radiative seesaw model with U(1)b-l gauge symmetry which is 
an extended model of the minimal HTM motivated from the type II seesaw mechanism 12]. 
The particle contents are shown in Table HT1 We add three right-handed neutrinos N c , vector- 
like SU (2) i doublet leptons L' and L' c , an SU (2) L triplet scalar A, an SU {2) L doublet scalar 
rj and B — L charged scalar \ to the SM, where r] and A do not have VEV. The Z 2 parity is 
also imposed so as to forbid the undesired terms. As a result, the neutrino mass is obtained 
not through the one-loop level (just like Ma- Model 13j) but through the two- loop level and 
the stability of DM candidates can be assured. Where we define V = (N4, E4) 44]. 

The renormalizable Lagrangians for Yuakawa sector and Higgs potential are given by 



•^Yukawa 
— -^Higgs 



^e c a L p + ytfN°L ,c + y 5 A L c s ir 2 AL' + yf X N c a N c b + ML'L' C + h.c, (II.l) 
2$t$ + mffirj + mlx ] X + m4Tr[A f A] - ^[$ T zr 2 A t ?7 + h.c] 
+ Ai($t$) 2 + \ 2 (rfv) 2 + H&$)(V ] V) + A 4 ($ t ^(r] t $) + A 5 /2[($ t r/) 2 + h.c] 
+ \{X ] X? + A 7 ( X t x)($ t $) + MX ] X){V ] V) + A 9 Det(A t A) 
+ Axo^rfA+A]) 2 + ai Tr[ti&\(rfri) + a 2 Ti[/^ r^) 

+ 6 1 Tr[A t A]($ t $) + 6 2 Tr[A t 7: i A]($ t r l $) + Cl Tr[A^ A]( X ^x), (H.2) 

where a, (3, 7, 5, a, b are the flavor indices. In the scalar potential, the couplings Ai, A2, Ae, 
and Ag have to be positive to stabilize the potential. Notice here that one straightforwardly 
probes the term &N C L, which derives neutrino mass at tree level, cannot be forbidden by 
any Abelian symmetries. If we introduce non-Abelian discrete symmetries 45|, |46[, one finds 



some successful groups such as Tj |47H50| and A (27) 5lN53| to forbid the term, remaining 
that 77 and A are one generation 3 . It implies that such an extension does not affect on 
analyses of dark matter and Higgs phenomenology focused on the inert scalar bosons. So 



3 We show in the appendix C how to realize. 
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we stay to analyze our model with flavor independent way hereafter for simplicity. 
The scalar fields x-, V an d A can be parameterized as 
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(II.3) 
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(II.4) 



where v is VEV of the doublet Higgs field $ satisfying v 2 = l/(y/2Gp) — (246 GeV) 2 , v' is 
that of the singlet Higgs field Xi an d subscript W denotes weak eigenstates. In Eq. (1II.3j) . 
G and G° (%/) are (is) the Nambu-Goldstone bosons (boson) which are (is) absorbed by 
the longitudinal component of W ± and Z (additional U(1)b-l gauge boson). The VEVs of 
A and r\ are taken to be zero because of the assumption of the unbroken Z 2 symmetry. 

In this model, the Z 2 -even fields $ and \ cannot be mixed with the Z 2 -odd fields A and 77, 
so that the mass matrices for the component scalar fields from $ and x an d those from A and 
77 can be separately considered. By inserting the tadpole conditions; m\ = — \±v 2 — Xjv' 2 /2 
and 7773 = — Xqv' 2 — Xjv 2 /2 (which is exactly the same as the results in Ref. |28|), the mass 
matrix for the CP-even scalar bosons in the basis of (0°, x°) is calculated as 



M\<P°,x°) 



2\ lV 2 X 7 vv' 
X 7 vv' 2\ G v' 2 



(11.5) 



The mass matrices of the ^2-odd scalar bosons are respectively obtained in the basis of (A^, 
77^), (Re A , ReT/°) and (ImA , Imr/°) by 



A4 2 (ImA°,Im77°) 
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M 2 (ReA° J Rer ] " 



777 (ReA°) fiv/V2 
fiv/V2 m 2 (Re7/°) 

(II.6) 
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where we define as follows: 

m 2 (A±) = m 2 + l Cl v' 2 + h lV 2 , (11.7) 

m 2 (ReA°) = m 2 (ImA°) = m 2 (A± ) + -v 2 , (11.8) 

™ 2 (Vw) = ™l + + ^ V/ 2 , (H.9) 

m 2 (Re V °) = m\ + ^\ 8 v' 2 + ^(A 3 + A 4 + \ 5 )v 2 , (11.10) 

m 2 (lm V °) = m 2 2 + ^X 8 v' 2 + ^(A 3 + A 4 - X 5 )v 2 . (11.11) 

The mass eigenstates for the Z 2 -even Higgs bosons and those for the Z 2 -odd scalar bosons 
can be defined by introducing the mixing angles a, /3, 7 and 5 as 



R((5) 





, cos 9 — sin 9 , 

with R(9) =| I , (11.12) 

sin 9 cos 9 

where h can be regarded as the SM-like Higgs boson. These mixing angles are expressed as 

t 2 A7W ^ 2^ f J,v 

Ai?; 2 -A 6 ?/ 2 ' m 2 (A^) -m 2 (?7^)' 

^ m 2 (ReA°) — m 2 (Rer/°) ' m 2 (ImA°) — m 2 (Im?7 ) ^ 

The mass eigenvalues are calculated by using the mixing angles given in Eq. (III. 13[) and the 
mass matrices given in Eqs. ( 111.51) and (III. 61) as 

R T (a)M 2 (^, X °)R(a) = \ ^ ° V R T ([3)M 2 (A± , r^)i?(/3) = ( ^ + ° 

m 2 H j y m 2 + 

R T (l)M 2 (A% v ° R )R( 1 )=[ ml « ° ) , R T (6)M 2 (AW I )R(5)= I ™ l < ° 

V </ V < 

(11.14) 

The doubly-charged triplet scalar bosons do not mix with others, and those masses are given 
by 

m 2 A++ =m 2 (A±)-^b 2 v 2 . (11.15) 
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FIG. 1: Neutrino mass generation via two-loop radiative seesaw. The particles indicated by a red 
font have the opposite %2 charge to those by a black font. 



We note that there is a characteristic relations for the mass spectrum among the triplet- 
like scalar bosons A ±=t , A^, A° R and A® in the limit of \x — > as 

m 2 A ++ - m\+ = m\+ - m 2 A o, (11.16) 
m 2 A0 = m 2 AR = m 2 Ai . (H.17) 

Therefore, two of four mass parameters for the triplet-like scalar bosons are determined by 
using above two equations. The same relations also appear in HTM when the lepton number 



violating coupling constant is taken to be zero [54]. 



B. Neutrino mass matrix 

The active neutrino mass matrix (depicted in Fig[l]) through two loop contribution is 
nven by 

(M] (yAUy„)l(M Nk )M 2 ( yA )p r d* P r A l i 

1 v>a P * I (27r) 4 7 (2tt) 4 (p 2 - M 2 ) 2 (p + q) 2 - M 2 Nk 



sin 2 7 cos 2 7 sin 2 5 cos 2 5 



p 2 — m 2 R p 2 — rn 2 AR p 2 — m 2 p 2 — m 2 A[ 
cos 2 7 sin 2 7 cos 2 5 sin 2 5 

X ~2 - + ~2 2 2 2 2 > ( IL18 ) 

q — rn nR Q — m A R 1 ~ m m Q ~ m A/ 
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FIG. 2: Lepton flavor violation. 



where 



d 4 p f d 4 q 



(2tt) 4 J (2tt) 4 (p 2 - M 2 ) 2 (p + g) 2 - p 2 - m 2 g 2 - m 2 

(3^i d1 / V- *>'(i- (IU9) 

M 2 H 

z z — z z z — z z — 1 



g 2 = -= M 2 + — m 2 + -m 2 . (11.20) 



As can be seen in the above equation, we can reproduce observed neutrino masses ~ 0(0.1) 
eV and the mixing data because of many parameters. 



C. Lepton Flavor Violation 



We investigate the LFV process £ a — > £pj (£ a , 1$ = e,fi,r) as shown in Fig. [2j The 



experimental upper bounds of the branching ratios are B — > e^) < 2.4 x 10~ 12 55], 
B(t -»■ fi'y) < 4.4 x 10" 8 and B (r -> e7) < 3.3 x 10~ 8 js6 1 . The branching ratios of the 
processes £ a — > ^7 are calculated as 



B(£ a -> ^7) 



3a e 
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{£ a ^£pi7pv a) , 
(11.21) 

where a cm = 1/137, B (ji —> eV^i/„) = 1.0, B (r — > eV^v T ) = 0.178, B (r — > [ivZv T ) = 0.174, 



Gf is the Fermi constant and the loop function ^(x) is given by 

1 — Qx + 3x 2 + 2x 3 — 6a; 2 In x 



Fox) 



6(1 -x) 



(11.22) 
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DM 
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FIG. 3: Dominant annihilation channel of DM. The left figure is for the fermionic DM; N c , and 
the right one is for the bosonic ones; rjR and A#. 

The fi — y process gives the most stringent constraint. The mixing angle sin/3 directly 
does not contribute to the neutrino mass, however, sin 2(3 should not be zero to retain fiv ^ 
in Eq. flILT3|) . 

III. DARK MATTERS 

We discuss DM candidates in this section. We have six DM candidates in general; that is, 
?7r(7), A#m, the lightest one of N c , and A^ 4 . However N4 cannot be the candidate because 
annihilation via Z boson gives too large cross section to obtain the observed relic 
density Vth 2 ~ 0.11 as well as too large scattering cross section in the direct detection 
search [?| 4 . Moreover, we restrict ourselves that only the real part of and are 
considered as a DM candidates, since the DM property of the imaginary part is more or less 
the same as the real one. Hence we analyze three DM candidates; 77??, A#, the lightest one 
of iV c , below. Hereafter we symbolize the DM mass as ttidm- 



In Ref. [4^, l57| , iV 4 is considered as a DM. 
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A. Fermionic Dark Matter 



We discuss a fermionic DM candidate N°, assuming the following mass hierarchy Mi < 
M2 < M3 for the right-handed neutrinos N?. Notice here that the DM mass range be less 
than M = 100 GeV to avoid the too short lifetime of DM. 

WMAP : At first, we analyze the DM relic density from WMAP. We have only the s- 
channel process via the Higgs bosons h/H as shown in the upper panel of Fig. |3j The 
effective cross section to / / /W + W~ /2Z° is given as 
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tDM 1 X ~ ^ " "** V2 ~ ^ " ^ 1 Ui > v ' (IIL1) 

where V = W ± ,Z° and we neglected the light quarks contributions and put only bottom 
quark contribution; that is, / ~ b, and i> re i; v Ie \ ~ 0(0.2) 58| . is the relative velocity of 
incoming DM. Here we neglect the contribution of H, assuming the mass is enough heavy. 
The SM-like Higgs mass is fixed to nih = 125 GeV. Notice here that the channel of VV 
is opened only if my < ttidm- in m^ < m^/2 < muM, the total decay width of h is 
= 4.1 x 10" 3 GeV [59(. Moreover, in rooM < mh/2, the channel h — > 2DM is also added 
and given by 



T h (h -»■ 2DM) 



irih f m-DM sin a 
8tt V v 1 



1 



/ m DM \ 
\ m h ) 



3/2 



(III.2) 



which is known as an invisible decay and recently reported by the LHC experiment that the 



60] 



branching ratio B inv is excluded to the region 0.4 < By 

Direct Detections : Let us move on to the discussion of direct detections. Our DM interacts 
with quarks via Higgs exchange. Thus it is possible to explore DM in direct detection 
experiments like XENON100 [7J. The Spin Independent (SI) elastic cross section ogi with 
nucleon N is given by 

Mdm / (Vs) u m N sma cos a 



m 2 n 



V2i 



f 



q 



(111.3) 
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where /x DM = ( 



■in 



DM 



+ m 



N ls j 1 is the DM-nucleon reduced mass and the heavy Higgs con- 



tribution is neglected. The parameters which imply the contribution of each quark to 



nucleon mass are calculated by the lattice simulation |6l|, |62| as 



fl = 0.023, f p d = 0.032, ff = 0.020, 

f: = o.oi7, ft = o.o4i, /; = 0.020, 



(III.4) 
(III.5) 



for the light quarks and fg = 2/27 ^1 — ^2 q<3 f^ for the heavy quarks Q where q < 3 
implies the summation of the light quarks. The recent another calculation is performed in 
Ref. 



Considering all the above constraints, we find that the allowed region is sharp at around 
m h /2. 



B. Bosonic Dark Matters 



We discuss bosonic DM candidates t}r and Ar 64|. Notice here that the DM mass range 



be less than M\y = 80 GeV to satisfy the constraint of the antiproton no excess reported by 



PAMELA 



65| as well as WMAP 



WMAP : At first, we analyze the DMs relic density from WMAP. We have two annihilation 
modes; t and u channel of 2r)n/An — > N^E^) — y 2v(££) and s-channel of 2r)n/ Ar — y h/H — y 
ff 6 . However since the t and u channel does not have s wave contribution, the dominant 
cross section is given only through the s channel 3l| as shown in the lower diagram of Fig. |3j 
The effective cross section to // is given as 
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(III.6) 

(III.7) 
(III.8) 



5 If the charged vector boson channel are open, the cross section is too large to satisfy the observed relic 
abundance. 

6 We neglect the H contribution. 
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When our DMs are less than m^/2, the invisible decay width is give by 31] 



T v R /Anf h 2DM) ~ i-* '-J I - Am 2 /m 2 h . (III.9) 

167rm h V ^ 1Vi 

Direct Detections : In the direct detection, the spin independent elastic cross section agi 
with nucleon N is given by 

2 / \ 2 



(T 



WA fl = Mdm m ^ [ A^ /Afl cos a 



where //dm and has been defined in the fermionic part. We here comment on the scenario 
based on the neutral component of the inert triplet field being the DM candidate, which has 
been discussed in Ref. [63]. Such a scenario is severely constrained by the direct detection 
experiments because of the Z boson exchanging contribution if the CP-even scalar boson 
and the CP-odd scalar boson from the triplet field are degenerate in mass. According to Ref. 
[63], the mass of the DM candidate has to be around 2.8 TeV in the model with Y=l inert 
triplet field in order to satisfy both the WMAP data and the direct search data. However, 
this result cannot be applied to our model, because we can take a mass splitting between 
A° R and A 5 due to the mixing between 77 and A. Thus, we can avoid the Z boson exchanging 
contribution to the direct search experiments. 

Considering all the above constraints, we find that the allowed region is the same as the 
fermionic case; that is, around m^/2. As a result, we assume to analyze that DM mass 
(especially A ) is in the next section. 



IV. HIGGS PHENOMENOLOGY 
A. Electroweak precision observables 

We discuss the constraints of the Higgs parameters from the electroweak precision ob- 



servables; i.e., the Peskin-Takeuchi S, T and U parameters [66j. If the mixing angle a is 
taken to be zero, the new physics contributions to the S and T parameters can be separated 
from those from the SM, so that we here assume a = for simplicity. Then the new physics 
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contributions to S and T are calculated as 

'4 Q (m|) - 4°(0) E| 3 (m|) - E 33 (0) 



'S'tirw — 167T 



(IV.l) 



T new = [£ 33 (0) - 5#(0)] , (IV.2) 

where S^ 1 , S 3 ^ and S 3 ^ 3 are obtained by calculating the 1PI diagrams of the gauge boson two 
point functions at the one-loop level whose analytic expressions are given in Appendix A. The 
experimental bound for the deviations in the S and T parameters from the SM prediction 
with the Higgs boson mass to be 117 GeV are given by fixing AU = as 

AS = 0.04 ± 0.09, AT = 0.07 ± 0.08, (IV.3) 

with the correlation factor to be 88% 



561 ] . The prediction in our model for AS" and AT are 
calculated by S ncw - (Ssm - S su \m h =ii7 Gev) and T ncw - (T SM -T SM |m h =ii7 Gev), respectively, 
where Ssm (Ssu\m h =in Gev) is the prediction for S in the SM with the SM Higgs boson mass 
to be m h (117 GeV), and T SM (Tg M | mh=117 GcV ) is those corresponding T parameter. We 
note that the exotic lepton V does not contribute to the S and T parameters, because both 
the masses of the component lepton fields (N± and E4) are given by M, and this field is 
introduced as the vector-like, so that both the custodial symmetry and the chiral symmetry 
are not broken by this field. 

In Fig. HI S-T plot is shown in the case where = m VR = m Vl = 300 GeV, and all the 
mixing angles are taken to be zero, so that the masses of the scalar bosons from the triplet 
field are determined by fixing two parameters; i.e., m&o and Am(= — m^o). Each 
dotted curve shows the prediction of the S and T parameters in the cases with m^o = 63 
GeV (m/j/2), 100 GeV and 200 GeV. The interval of each dot indicates the increment of 
Am{= m A + — m A o) to be 3 GeV from right to left. The inside (outside) ellipse indicates 
allowed region with the 68% (95%) confidence level by the electroweak precision data. It 
is seen that that the T parameter is getting larger values when Am is taken to be larger 
values. When m A o is taken to be 63 GeV, 100 GeV and 200 GeV, the allowed maximum 
value for Am with 95% confidence level is about 15 GeV, 23 GeV and 30 GeV, respectively. 



11 




"°'3).2 -0.1 0.1 0.2 

s 

FIG. 4: Prediction of the S and T parameters for the case with m n + = m rj o(= m rjR = m Vl ) = 300 
GeV. All the mixing angles (a, /3, 7 and 5) are taken to be zero in this plot. Each dotted curve 
is shown the results in the cases with m^o to be 63 GeV, 100 GeV and 200 GeV, where the 
interval of each dot indicates the increment of Am(= — m^o) to be 3 GeV from right to 

left. The regions within the blue (red) ellipse are allowed at the 68% (95%) confidence level by the 
electroweak precision data. 



B. Signature of A ±=t at LHC 

We discuss how our model can be tested at collider experiments. We focus on the signa- 
ture from the doubly-charged scalar bosons A ±=t , because appearance of A ±=t is one of the 
striking properties of the model. In order to focus on the phenomenology of the triplet-like 
scalar bosons, we assume that the mixing between A and 77 are taken to be quite small 
(/3 ~ 7 ~ 5 ~ 0), and the masses of rf 1 and 77° are much heavier than those of the triplet-like 
scalar bosons. 

There are an indirect way and a direct way to identify existence of A ±=t . The farmer way 
is measuring the deviations in the event rate for the Higgs boson decay channels from the 
SM values. In particular, the Higgs to diphoton mode h — > 77 is one of the most important 
channels for the SM Higgs boson search at the LHC because of the clear signature. The 
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current signal strength for this channel is 1.65±0.24 at the ATLAS l|,|2| and 1.6±0.4 at the 
CMS The decay rate of h — > 77 can be modified by the loop effect of A ±=b . Contributions 



from doubly- charged scalar 
several papers in the HTM 



Dosons to the h — > 77 mode have already been analyzed in the 



67 



68[. The signal strength for h — > 77 can be larger than 1.6 



in the case with m^±± to be smaller than about 200 GeV 68[ without contradiction with 



the constraints from the vacuum stability [69| and the perturbative unitarity [69|, |70|. The 
prediction of the deviation in the decay rate of h — > 77 in our model is almost the same as 
that in the HTM as long as the contributions from rf 1 are neglected. 

Discovery of A ±=t can be a direct evidence of our model. We consider the case where 
a lighter neutral scalar boson from the triplet field (Ar or Aj) is assumed to be DM. In 
that case, the neutral scalar boson should be the lightest of all the triplet-like scalar bosons; 
i.e., m^o < m&+ < m&++ to guarantee the stability of DM, and its mass is around the half 
of the Higgs boson mass to satisfy WMAP data and direct detection experiments. In 
addition, the mass differences among the triplet-like scalar bosons are restricted from the S 
and T parameter as shown in Fig. HJ e.g., Am is constrained to be less than 15 GeV in the 
case with m&o = 63 GeV. This implies that the upper limit for m^++ is about 90.5 GeV by 
using the mass relation given in Eq. (III.16j) . 



A search for doubly-charged Higgs bosons has been done at LEP [7l|, Tevatron |72[ and 



LHC 



73]. All these searches have been done under the assumption where doubly- charged 



Higgs bosons decay into the same sign dilepton. The most stringent lower bound for the mass 
of doubly-charged Higgs bosons is about 400 GeV given at LHC [t3 | . In our model, A ±=t 
cannot decay into the same sign dilepton associated without any other particles, because 
they are the Z 2 odd particles. In that case, the mass bound given at LHC cannot be applied 
to that in our model. 

First, we discuss the decay of A ±=t in the scenario based on the lighter neutral component 
of A assumed to be DM. Basically, there are two decay modes of A ±=t : (1) the same sign 
dilepton decay through the Yukawa coupling y& (A ±=t — > Ef^) 7 with to be e ± , ^ or r*, 
and (2) W boson associated decay through the gauge coupling constant (A ±=t — > A ± W ± *). 
The decay branching fractions for (1) and (2) are determined by the magnitude of y&, Am 



7 The magnitudes of the branching fractions of A ±=t in the same sign dilepton modes: Efe^, and 
E^t^ depend on the value of y q A . In the following discussion for the collider phenomenology, we do not 
specify the flavor of £ ± . 
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M [GeV] 



FIG. 5: The decay branching ratio of A ±=l: as a function of the exotic lepton mass M for the case 
of m A ++ = 90.5 GeV, m A + = 78 GeV. The solid and dashed curves respectively show the cases 
with y A = 0.1 and 0.01. 



and the mass of the exotic lepton M. The formulae for the decay rates for these channels 
are given in Appendix B. 

In Fig. [5j the branching fraction of A ±=t is shown as a function of M. We take m A ++ = 
90.5 GeV and m A + = 78 GeV which correspond to the case with m A o = 63 GeV and 
Am = 15 GeV. The Yukawa coupling y A is taken to be 0.1 and 0.01. It is seen that the 
main decay mode is changed from A ++ — > E\i + to A ++ — > A + W + * when M is getting 
larger values. For example, 50% of i3(A ++ — > E^£ + ) can be obtained in the case of M ~ 89 
GeV (84 GeV) and y A = 0.1 (0.01). 

In the following, we discuss the case where A ±=t mainly decay into the same sign dilepton 
(A±± £ ± Ef). In that case, the exotic lepton mass M should be between m&++ and 7«a+, 
otherwise there is no possible decay channel for Ef. As an benchmark scenario, we take the 
following mass spectrum and the coupling constant which is allowed from the electroweak 
precision data and also the LFV data: 

m A++ = 90.5 GeV, m A+ = 78 GeV, m A o = 63 GeV, M = 85 GeV, 

y A = 0.1. (IV.4) 

In this parameter set, the branching fraction of A ±=t — > Ef^ is about 99%, and that of 
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— > W /± A° and Ef — y A ± u are 100%, where A is A^ or Aj. Thus, the decay process of 
A ++ is expected as follows: 

a ++ -> e + E+ ->■ £+a+z/ -> £ + iy+A°z/ ->■ e + e + A°uu, (iv.5) 

so that the final state contains the same sign dilepton and the missing energy. 

We then discuss the signal of A ±=t at LHC in the parameter set given in Eq. (IIV.4I) . At 
LHC, A ±=t can be mainly produced via the Drell-Yan processes: qq — > ^y*/Z* — y A ++ A _ " 
and qq' — > W^* — > A ±± A =F . Thus, the signal events are expected to be 8 

qq -> A++A" -»• (£ + W + *vA°)(rW-*vA°) -> ffrt| T , (IV.6) 
gg' -> A ±± A =F -> (£ ± W ± VA°)(^ T *A°) £ ± £ ± £^^ T . (IV. 7) 

The cross sections of cr(pp — > A ++ A ), cr(pp — > A ++ A~) and a(pp — > A A + ) are 
respectively evaluated as 566 fb, 889 fb and 494 fb with the collision energy to be 8 TeV by 
using CalcHEP jlJ] and CTEQ6L parton distribution functions. The cross sections for the final 
states expressed in Eqs. (1IV.6j) and (1IV.7|) are obtained as 63 fb and 154 fb, respectively. 



The same events can happen in HTM. The doubly-charged Higgs bosons if ±=t in HTM 
can decay into the same sign diboson if ±=t — y W ± W ± which is realized in the parameter 
regions where VEV of the triplet field va is larger than about 10~ 4 GeV, and if ±=t is the 
lightest of all the Higgs bosons from the triplet field. At the same time, singly-charged Higgs 
bosons H in HTM can decay into W ± Z in this case as long as the mass difference between 
H and H ± is not too large. Therefore, the same events as expressed in Eqs. (IIV.6I) and 
(1IV.7I) can appear in the following way: 

pp -»■ H ++ H— -»■ W + W + W~W~ £+£+£-£- # T , (IV.8) 
pp -»■ -»■ W ± W ± W T Z -»• PPFfyr, (IV.9) 

Furthermore, in the case where if ±=t are the heaviest among the triplet like Higgs bosons 
such like our scenario, the cascade decay of if ±=t can be dominant; i.e., H ++ — y W ± H ± —y 
W ± W ± H° (H° is a neutral component field of the triplet Higgs). In addition, v& < 10 -4 
GeV, H° mainly decays into neutrinos. In this case, the final states of the signal event can 



The heavier A ; e.g., An can further decay into Z* and Aj which is corresponding to DM. 
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FIG. 6: Invariant mass distribution for the £ + £ + system from the event pp — > £ + £ + £~£~^t in our 
model (red solid curve) and that in HTM with the diboson decay (black dotted curve) of and 
the cascade decay of (blue dashed curve). We take the mass of the doubly-charged scalar 

bosons to be 90.5 GeV and the collision energy to be 8 TeV. The solid, dashed and dotted curves 
are respectively shown the distributions for the dilepton decay in our model, cascade decay in the 
HTM and diboson decay in HTM. Number of events are assumed to be 10 4 in each distribution. 



also be 

pp -> H ++ H-~ -> (W + W + H°)(W-W-H°) -> £+£+£-£- # T , (IV. 10) 

pp -> H ±:k H T -> (W ± W ± H°)(W T H°) -> tl±F$ T . (IV. 11) 

The invariant mass distributions in the system of the same sign dilepton Me+e+ can be 
useful to discriminate between our model and HTM. In Fig. [61 the distribution for M(+(+ 
in the £ + £ + £~£~ Ifir system is shown in which the number of event is assumed to be 10 4 and 
collision energy to be 8 TeV. The mass of the doubly-charged scalar bosons in our model 
and also in HTM are taken to be 90.5 GeV. The red solid curve and the black dotted (blue 
dashed) curve are respectively represent the M1+1+ distribution from our model and HTM 
with the diboson decay (cascade decay) of H . It can be seen that the event number 
from the diboson decay of if ±=t is distributed in the wide region of M £ + £ +, and it takes 
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the maximum value at around the half of the mass of i7 ±=t ; i.e., M^+^+ ~ 45 GeV. On the 
other hand, the other two shapes of the distributions look like each other in which the event 
number is distributed in the small region of Mg+g+. However, the cross sections for the final 
states are different between these two events. In our model, one of the two leptons with the 
same sign comes from the decay of A ±=t , and the other one comes from the leptonic decay 
of the W boson, so that the cross section is calculated as a(pp — >■ A ++ A ) x B(W —> £v) 2 ■ 
On the other hand, in HTM with the cascade decay of if ±=t , all the leptons in the final state 
are obtained from the leptonic decay of the W boson, so that the cross section is evaluated 
by cr(pp — > H ++ H ) x B(W — > iu) 4 . The original cross sections of a(pp — > A ++ A ) and 
<j(pp —> H ++ H ) are the same each other as long as we take the masses of A ±=t and if ±=t 
to be the same. Therefore, the event number in HTM with the cascade decay scenario is 
smaller than that in our model by the factor of B(W —> lv) 2 = 1/9. 

We then conclude that the four lepton events with the missing energy £ + £ + £~£~$t from 
our model and HTM with the diboson decay and the cascade decay of i? ±=t may be able to 
be discriminated by using the M1+1+ distribution and the number of event. 



V. CONCLUSIONS 



We have constructed a two-loop radiative seesaw model that provides neutrino masses in 
a TeV scale theory. We have also studied DM properties, in which our model has fermionic 
(N c ) or bosonic DM (77°, A ) candidate with the same mass scale, which is at around 
from the constraint of WMAP and the direct detection search in XENON100. We have also 
discussed Higgs phenomenology at LHC, in which the neutral scalar boson from the triplet 
field is assumed to be DM. In that case, the mass of the doubly-charged scalar bosons A ±=t 
is constrained to be smaller than about 90.5 GeV by the electroweak precision data. In 
this scenario, A ±=t can mainly decay into the same sign dilepton with the missing energy 
where the exotic leptons E± appear the intermediate state of the decay process of A ±=t . We 
then considered how we can discriminate the signature of A ±=t in our model from that of 
doubly-charged Higgs boson if ±=t from HTM. We have found that the distribution of the 
same sign dilepton system and the magnitude of the cross section for £ + £ + £~£~I£t can be 
useful to distinguish between our model and HTM. 
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We list the analytic formulae for the 1PI diagram contributions to the gauge boson 
two point functions at the one loop level which are necessary to calculate the S and T 
parameters. We ignore the mixing between the Z 2 -even scalar bosons; i.e., a = 0. In that 
case, new physics contributions to the gauge boson two point functions can be separated 
from the SM ones. The new physics contributions are calculated by 



Appendix A: Gauge boson two point functions 



Y% w {p 2 ) 




+ Ac 2 B 5 (p 2 ,m A++1 m A +) + As 2 B 5 (p 2 , m A ++, m v +) 

+ (V2cpOy + sps y ) 2 B 5 (p 2 ,m A +,m AR ) + (V2cpc s + s f3 ss) 2 B 5 (p 2 ,m A +,m Al ) 
+ (c/3C 7 + A/2s /3 s 7 ) 2 J B 5 (p 2 ,m r? +,m^) + (cpc s + V2spS S ) 2 B 5 (p 2 ,m v + ,m Vl ) 
+ (-V2cps y + s /3 c 7 ) 2 J B 5 (p 2 ,m A +,m^) + (->/2cps s + s f3 c 8 ) 2 B 5 (p 2 ,m A +,m Vl ) 
+ (-a/2s/3C 7 + cps 1 ) 2 B b (p 2 ,m n +,m AR ) + (-y/2spc s + cpssf 'B b (p 2 \m n + ,m Al ) , 



(A.l) 



= i6^K jB3(p2,M ' M) 

+ 4B 5 (p 2 , m A ++ , m A ++ ) + B 5 (p 2 ,m A +, m A + ) + B 5 (p 2 , m v + , m v + ) , 



(A.2) 



Zp(p 2 ) 



+ 2(1 - 2s 2 v )B 5 (p 2 ,m A ++,m A ++) + -(sj - 2s 2 w )B 5 (p 2 , m A +, m A +) 
I 

+ ~ 2s w) B s(p 2 i m r,+ i m r)+) , 



(A.3) 
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J_gJ 

16tt 2 4 



32p 2 (^-s 2 w + s 4 w )B 3 (p 2 ,M,M) 



+ 4(1 - 25^) B 5 (p ,m A ++,m A ++) + (s^ - 2s w ) B 5 (p ,m A +,m A +) 
+ (cj -2s 2 v ) 2 B 5 (p 2 ,m ri +,m v +) + 2c 2 p s 2 ^B 5 (p 2 , m A +, m v +) 



{2c^c s + s^s s ) 2 B 5 (p 2 , m Afl , m A/ ) + (c 7 c<5 + 2s 1 s s ) 2 B 5 (p 2 , m VR ,m 



'ii , 



+ (-2c 7 s,5 + s y cs) 2 B 5 (p 2 , m Afl , m 



-2s 7 c,5 + c^) 2 ^^ 2 , m^, m A/ ) , (A.4) 



where eg = cos 6* and = sin#. In the above equations, B 3 (p 2 ,m 1 ,m 2 ) and B 5 (p 2 , m J _ L _ m 2 ) 
functions [3] are respectively expressed in terms of the Passarino-Veltman functions 76 1 by 



B 3 (p 2 ,mi,m 2 ) = -^(pV^m.) - B 21 (p 2 ,m u m 2 ) 
B 5 (p 2 , mi, m 2 ) = A(m 1 ) + A(m 2 ) - 4B 22 (p 2 , mi, m 2 ) 



(A.5) 
(A.6) 



The functions E^ , E 3 ^ and E^ 3 are given in terms of above the gauge boson two point 



functions by 



^E^ + E 



sw 



V 33 — 

# 2 



2 y,ZZ 



2sy/C\\/Yij, 



Appendix B: Decay rates 



(A.7) 



The decay rates for the doubly- charged scalar bosons A ±=t are calculated as 



r(A ±± -»• ^e ±n 



m; 



M 2 



m 



A+^ 



m 



A+^ 



A 1/2 



m; 



M 2 



m 



A+^ 



m 



A++ 



T(A ±± A ± W ± ) = cos 2 /3A 3 / 2 f 

16tt m^ W A ++ m A++J 



r(A ±± -> A ± W 



16n mly 
9g 4 m A + 
128vr 3 

2 ™3 



cos 2 /?G 



r 



A+ 



m 



T(A ±± ^W*) = -^!%±sin 2 /3A 3 / 2 ( ^$ 



m 



u 



16tt mly 
9g 4 m A ++ . 



\ m A++ m A++ 



m 



r(A ±± -+ ^w^*) = \ Z~* sin2 £ g . 



128tT 



(B.l) 
(B.2) 
(B.3) 

(B-4) 
(B.5) 
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where X(x, y) and G(x,y) are the phase space functions which are given by 

X(x, y) = 1 + x 2 + y — 2x — 2y — 2xy, 
1 



G(x,y) 



12y 



2 (-1 + x) 3 - 9 (-1 + x 2 ) y + 6 (-1 + x) y 2 

1 + x — y 



+ 6(1 + x - y)y\/-X(x,y) 
- 3 [1 + (x - yf - 2y] ylogx 



tan 1 



tan 



(B.6) 



— 1 + x + y 
y/-X(x,y) 



(B.7) 



Appendix C: Flavor symmetry 



Particle 


L 


e c 


N c 


V 


L' c 


A 




V 


(i) 

Xi 


x\ 


T 7 


3 


3 


3 


3 


3 


lo 




lo 


3 


3 



TABLE III: The particle assignments in T? symmetry. Here i runs 0-2. 

Here we show an example how to realize our model for the lepton sector by using non- 
Abelian discrete symmetry. The minimal extension is to introduce TV flavor symmetry 
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50|. Each of the field assignment is given in Table [LTT| where the other assignments are 
same. The extended Lagrangian to Eq. ( 111.11) is modified as 9 

3 333 



i=l 



i=l 



i=l 



i=l 



+ Vs (x?{N*N° + A^A 2 C ) + xf^iVT + N°N°) + (A^l + A 2 C A X C )) + h.c, (CJ 

where u = e 2t7T ^ 3 . Here notice that all the terms except yA c iV c are diagonal. It suggests that 
the observed neutrino mass and lepton mixing can be obtained only through the xA c A c 
term. After the B — L spontaneously breaking; {x!i) — 1$ an d {X2) = v 2 > the right- 



9 The chargcd-lcpton sector has to be improved in order to forbid the universal Yukawa coupling of the 
SM-likc Higgs boson, since it has been ruled out by the current Higgs boson search data at LHC. Straight- 
forward ways to solve it is to change the flavor symmetry group and/or introduce some additional Higgs 
fields. 
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handed neutrino mass matrix is given by 



/ (2) 



y a s 



(3) 



o 





V 



vl 



v[ 1] J 



( vi 2) \ 



y b s 



V 



(3) 



V 



(1) 



V 



(2) „« 







(C.2) 



/ 



As a result, we can easily find the observed neutrino mass difference and their mixings by 
controlling each of the VEV. In the Higgs potential, the modifications are as follows: 



2 3 

yt Y V V v t(j) y (i) 
i=l j=l 



(C.3) 
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